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ABSTRACT 

We report here on a sequence of 28 observations of the binary pulsar system 
PSR B1259— 63/SS2883 at four radio frequencies made with the Australia Telescope 
Compact Array around the time of the 2000 periastron passage. Observations made 
on 2000 Sep 1 show that the pulsar's apparent rotation measure (RM) reached a 
maximum of —14800 ± 1800 rad m~ 2 , some 700 times the value measured away from 
periastron, and is the largest astrophysical RM measured. This value, combined with 
the dispersion measure implies a magnetic field in the Be star's wind of 6 mG. We find 
that the light curve of the unpulsed emission is similar to that obtained during the 
1997 periastron but that differences in detail imply that the emission disc of the Be 
star is thicker and/or of higher density. The behaviour of the light curve at late times 
is best modelled by the adiabatic expansion of a synchrotron bubble formed in the 
pulsar/disc interaction. The expansion rate of the bubble ~12 km s _1 is surprisingly 
low but the derived magnetic field of 1.6 G close to that expected. 

Key words: radiation mechanisms: non-thermal — binaries: close — stars: emission- 
line, Be — stars: individual: SS2883 — pulsars: individual: PSR B1259— 63 — radio 
continuum: stars. 



1 INTRODUCTION 

PSR B1259— 63 was discov ered in 1990 in a surv ey of the 



Galactic plane for pulsars (Johnston et al. 1992a) and was 



subsequently found to be in a 3.4-y r, highly eccentric orbi t 
about a ~10 M Q Be star, SS2883 ([Johnston et al. 1992b|). 



It remains the only known radio pulsar with a Be star com- 
panion. The pulsar has a spin period of ~48 ms, a charac- 
teristic age of 0.33 Myr, and a moderate magnetic field of 

3.3 x 10 11 G. The dispersion measure (DM) of 146.7 cm 3 pc 
yields a distance of 4.5 kpc in the model of Taylor fc Corde s 
(1993), howe ver scintillation (|McClurc-Griffiths et al. 1997j ) 
and optical (Johnston et al. 1994) measurements indicate a 
distance closer to 1.5 kpc. The pulse profile, shown in Fig- 
ure El has two almost equal intensity peaks both of which are 
highly (and almost orthogonally) polarised, with an average 
interstellar rotation meas ure (RM) measured away from pe- 
riastron of +21 rad m~ 2 ( Manchester fc Johnston 1995 ). 

Optical observations ( Johnston et al. 1994| ) show that 
SS2883 is of spectral type B2e, and thus has a mass of ~10 
Mq and a radius R c ~6 Kq. Assuming a pulsar mass of 

1.4 Mq, the implied inclination angle of the binary orbit 
to the plane of the sky is 36°. The Ha emission line shows 
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Figure 1. Mean pulse profile at 1.4 GHz. Position angle is shown 
on top, and the total intensity (solid line), linear (dashed line) and 
circular (dash-dot line) polarisations are shown in the bottom 
panel. 



that the Be star's emission disc extends to at least 20 R c , 
just inside the pulsar's orbital radius at periastron. Timing 
measurements have shown that the disc of the Be star is 
l ikely to be highl y inclined with respect to the orbital plane 
( |Wex et al. 1998j ). 

The 1997 periastron was observed extensively by John- 
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ston et al. (1999, 2001) and a model for the transient emis- 
sion was proposed by Ball et al. (1999) following earlier work 
by Melatos et al. (1995). In brief, unpulsed emission began 
at T — 22 and lasted until at least T+ 100 (where T denotes 
the epoch of periastron) . The spectral index of the unpulsed 
emission was in the range —0.5 to —0.7, a value indicative 
of synchrotron radiation. Peaks in the light curve at T — 10 
and T + 20 coincide with the pulsar crossing through the 
emission disc of the Be star. The model therefore proposed 
that the synchrotron emission was generated in the shock 
between the relativistic pulsar wind and the outflowing disc 
of the Be star. Synchrotron losses, set by the magnetic field 
strength, then determined the decay time of the light curve. 
The decay time of the emission from the post-periastron 
disc crossing was substantially shorter than that of the pre- 
periastron decay. 

The model predicted a frequency dependence for the 
end of the transient phase and verification of this late-time 
behaviour was a motivation for the 2000 periastron observa- 
tions. Also, the improvements in the correlator and the data 
reduction software meant that the pulsar rotation measure 
(RM) and dispersion measure (DM) could be accurately ob- 
tained during this periastron passage. We therefore obtained 
data over 28 epochs from T - 46 to T +113 during the 2000 
periastron encounter. In section 2, we describe the observa- 
tions. In section 3 we detail the results obtained and discuss 
their implications in section 4. 



2 OBSERVATIONS 

Observations were made at the Australia Telescope Com- 
pact Array (ATCA). The ATCA is an east- west synthesis 
telescope located near Narrabri, NSW and consists of six 
22-m antennas on a 6 km track. On the ATCA, observa- 
tions can be made simultaneously at either 1.4 and 2.4 GHz 
or 4.8 and 8.4 GHz with a bandwidth of 128 MHz at each fre- 
quency subdivided into 32 spectral channels, and full Stokes 
parameters. The ATCA is also capable of splitting each cor- 
relator cycle into bins corresponding to different phases of a 
pulsar's period, and in our case the pulse period of ~48 ms 
was split into 16 bins. 

A total of 28 observations between T — 46 and T + 113 
were made using various ATCA configurations. Table |l] lists 
the relevant details. One of the ATCA primary flux calibra- 
tors, B0823-500 (flux densities of 5.5, 5.6, 3.1 and 1.4 Jy at 

I. 4, 2.5, 4.8 and 8.4 GHz respectively) or B1934-638 (14.9, 

II. 6, 5.8 and 2.8 Jy), was observed once during each ses- 
sion for ~5 min. Each pair of bands was switched every ~20 
min, with an observation of the secondary phase calibrator 
(B1251— 713; flux density ~1 Jy at all frequencies) for ~3 
min leading each integration. 

Data reducti on and analysis were carried out with the 
Miriad package (Sault & Killeen 1998) using standard tech- 



niques. After flagging bad data, the primary calibrator was 
used for flux density and bandpass calibration and the sec- 
ondary calibrator was used to solve for antenna gains, phases 
and polarisation leakage terms. The short baselines in the 
750C array typically were completely flagged out at the two 
lower frequencies, because of excessive interference. After 
calibration, the data consist of 13 independent frequency 
channels each 8 MHz wide for each of the 16 phase bins. 
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Table 1. Date, length of observation in hours, and array con- 
figuration for all observations of PSR B1259-63/SS2883 during 
the 2000 periastron passage. On Nov 7, only 1.4 and 2.4 GHz ob- 
servations were made. On Dec 8 time-binning operation was not 
functional. Periastron occurred on Oct 17.3. 



2.1 Total Flux Density 

To measure the total flux density from the system (i.e. the 
pulsar plus any unpulsed emission) we first collapsed the 16 
phase bins into one. Standard techniques of image inversion, 
cleaning and restoring were then performed. The flux density 
and associated errors were measured directly from the image 
with the routine imfit with the source assumed to be a point 
source at a fixed position given by the timing data (Wex 
1998|) . 



et al. 



2.2 Pulsar Flux Density 

The uv binned data were first corrected for dispersion 
with the routine psrfix. We assume a constant DM of 
146.7 cm _3 pc for all the observations. Although the DM 
is known to vary close to periastron, our time resolution 
of only 3 ms made the increased time delay largely irrele- 
vant. We then examined the pulse profile using and selected 
those bins which formed the off-pulse portion of the profile. 
The mean of the selected baseline bins is then subtracted 
from all the phase bins for every frequency channel and for 
each Stokes parameter using psrbl. This has the effect of 
removing all sources from the image apart from the pulsar. 
The flux density of the pulsar was then obtained directly 
from the visibility data (with uvflux) without the need 
for further imaging. This process is valid in the absence of 
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Table 2. Flux densities for pulsed emission and total emission (pulsed plus unpulsed) for all 29 observations with the ATCA. Flux 
densities and errors for continuum and pulsar are from IMFIT. Flux densities and errors for pulsar only are from uvflux. Pulsar flux 
densities are omitted when no significant pulses could be seen. 



sidelobes from any other source in the primary beam. The 
main source of error in this technique is determining which 
phase bins are off-pulse. This is especially true when the 
pulsar is weak and/or scattered. Generally we found that 
the error in the flux density given by the imaginary compo- 
nent in uvflux was equivalent to the value expected from 
thermal system noise considerations. Occasionally, a signifi- 
cantly larger measurement error was determined (indicated 
with bold type in Table ^) which was dominated by phase 
calibration errors. 



2.3 Dispersion Measure 

Even with our relatively crude time resolution of ~3 ms 
per bin we are able to measure the DM across the band 
at 1.4 GHz with reasonable accuracy. To achieve this we 
determined the centroid of each of the two pulse components 
for each of the 13 frequency channels. The time delay in 
seconds for a given channel, c, relative to the first channel 
is given by (see e.g. Manchester & Taylor 1977) 



on the DM is about one order o f magnitude higher t han 



1,, = 4149.94 DM | \ - X 



(1) 



where v is the frequency of the channel in MHz. The DM 
can therefore be obtained by a straight-line fit to a plot of 
dispersion delay versus wavelength squared. The error bar 



from 'conventional' pulsar timing (Johnston et al. 2001) 



2.4 Circular Polarisation 

As the natural modes of the wind of the Be star and the in- 
terstellar plasma are circular, circularly polarised flux is not 
affected by propagation. We therefore used the measurement 
of the circular polarisation under each pulse component in- 
dependently to verify our polarisation calibration. As before, 
the flux in Stokes V for two separate phase bin ranges corre- 
sponding to the two pulsar components was obtained from 
the visibility data via uvflux. 



2.5 Linear Polarisation and RM 

UVFLUX was used to measure Stokes Q and U in each chan- 
nel for each phase bin and/or combinations of channels and 
phase bins. In practice, in order to boost signal to noise, we 
combine the phase bins together which correspond to each 
pulse component as Stokes Q and U do not vary rapidly 
through the pulse (see Fig 1). The position angle, cj>\, of 
each pulse component at a given observing wavelength, A, 
can be then be computed through 



(/)x = 0.5 x tan 1 {Ux/Qx). 



(2) 
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Figure 2. Light curves for the unpulsed emission (top) and the 
pulsar (bottom) from the 2000 periastron data. For clarity the 
error bars have been omitted. 



The error on the position angle is given by 



A0 = 0.5 £ /a/Q 2 + U 2 



(3) 



when both the Q and U measurements have the same rms 
error e. Note that as the linear polarisation has a positive 
definite bias, this must be subtracted before the error is 
computed. The position angle is related to the RM by 



= RM A 



(4) 



Hence, plotting position angle against A 2 should yield a 
straight line with slope RM. This procedure was carried out 
for each of the two components independently. 
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Figure 3. Spectral indices for the unpulsed emission. 



days later. At T — 18 the evidence for pulsed emission at 
any frequency is marginal although it cannot be ruled out. 
Following this date there is no further detection of pulsed 
emission until T + 18. Post-periastron, there is a clear de- 
tection at the two higher frequencies at T + 18 and shortly 
after the end of the ATCA observation the pulsar was also 
detected at Parkes at 1.3 GHz. Following this, the pulsar 
was detectable again at all frequencies. 

It seems likely therefore that scattering dominates the 
detectability of the pulsar for ~20 days leading up to T — 
16. From this date until T + 18 the pulsar is behind the 
emission disc with respect to the observer and the pulsed 
emission is quenched by the high optical depth. The pulsar 
then appears to emerge from the disc rather abruptly; the 
scattering is marginal although optical depth effects are still 
visible at 1.4 GHz until T + 25. Unfortunately, the poor 
time resolution of our data precludes measurement of the 
scattering parameters. 

In 1997 the pulsar was detected with the ATCA until 
T — 21 at all frequencies and at Parkes until T — 18. Post- 
periastron it was detected from T + 16 onwards. The effects 
of the wind on the pulsed emission therefore occurred earlier 
and lasted longer in 2000 than in 1997. 



3 RESULTS 

The flux densities and their associated errors from the ob- 
servations listed in Table jj] are given in Table |^. The table 
first lists the date of the observation, followed by the time 
from periastron in days. The rest of the table contains the 
flux density measurements of the pulsar (when detectable) 
and the total flux density (pulsed plus unpulsed) obtained at 
each of the four frequencies using the techniques described 
in the previous section. 



3.1 Pulsed Emission 

The pulsar fluxes are shown in the lower panel of Figure ^| 
Before periastron, the pulsar was detected at all frequencies 
at T — 46. At T — 32 the pulsar was scattered at 2.5 GHz and 
not detected at 1.4 GHz, at T — 24 scattering was evident 
at 4.8 GHz and the pulsar not detected at the two lower 
frequencies. However, at T — 22 the pulsar was again de- 
tected at 2.5 GHz but was scattered at that frequency two 



3.2 Unpulsed emission 

The top panel of Figure |^ shows the light curve for the 
unpulsed emission. Unpulsed emission appears to be present 
at 1.4 GHz during the first observation at ~ T — 46 and 
remains roughly constant until ~ T — 20. The flux density 
then rapidly increases (to 40 mjy at 1.4 GHz) towards the 
first peak at ~ T — 6. A decrease with a similar slope follows, 
between ~ T — 6 and T + 2. The flux density plateaus (at 
~ 21 mjy at 1.4 GHz) at all frequencies between T + 1 
and T + 13, and then rises again to a peak at ~ T + 22. 
The rise to the second peak is slower than the first, and the 
subsequent decay is rather shallow. T + 70 shows a sudden 
drop in the flux density, subsequently however the unpulsed 
emission remains clearly present at the two lower frequencies 
during our final observation at T + 113. This is similar to 
data taken in 1994 which shows unpulsed emission p resent 
at 0.84 GHz until at least T +132 ( [Johnston et al. 1999| ). The 
data do show, however, that unpulsed emission is not present 
at 8.4 GHz from T + 63 onwards and has also vanished at 
4.8 GHz by T+ 113. 
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Figure 4. Position angle versus A 2 across the 4.8 GHz band on 
September 1 for the first (top) and second (bottom) components. 
The position angle separation between the components is near 
flO°, as expected. The straight lines show the fit to the data 
from which the RM can be derived. 



3.3 Spectral indices 

The flux density at a frequency v is given by S v = Cv a 
where a denotes the spectral index. We measure the spec- 
tral index for the non-pulsed emission using flux densities 
obtained at each of our four observing frequencies. We apply 
a simple linear fit to log S - log v space and assume sym- 
metrical errors in log S. The spectral indices are shown in 
Figure ^ and are consistent with a value of ~ —0.7 through- 
out. There is no evidence for optical depths effects, at least 
at freque ncies above 1 GHz. T his is similar to the behaviour 
in 1997 jjohnston et al. 199S| ). 

We do not attempt to measure the spectral index of 
the pulsed emission. At 8.4 GHz, diffractive scintillation sig- 
nificantly affects the observed flux density. The sc intillation 
bandwidth of ~70 MHz and timescale of ~10 min ( McClure- 



Griffiths et al. 1997) match well with our observing band- 
width of ~120 MHz and time of 20 min. Furthermore, there 
is clear evidence for optical depth effects at low frequency 
post-periastron. In particular, the flux density at 1.4 GHz 
appears to increase steadily from T + 21 to T + 70 whereas 
this is not the case at the two higher frequencies. 



3.4 Polarisation and Dispersion Measure 

The circular polarisation was obtained for each pulse com- 
ponent for each observation, with errors of order ~ 5% at 
the lower frequencies, and ~ 10% at the higher frequencies 
due to the higher system temperature and lower flux density 
at these frequencies. The data are consistent with a mean 
value in the range —6% to —10% for the first component, and 
between —28% and —30% for the second component for all 
frequencies, in agreement with observations made at Parkes 
(see also Fig []]) . We are thus confident that our polarisation 
calibration and analysis techniques are correct. 

We attempted to obtain RMs and DMs for all days when 
the pulsar was not eclipsed using the techniques described in 
Section 2. These are presented in Table ^. The RM values are 
weighted averages of each of the two pulse components. In all 
cases the RMs agreed within the errors. Note that in general 
for RMs <2000 rad m -2 , only the two lower frequencies will 



show a significant position angle change across the bandpass 
whereas for RMs > 5000 rad m -2 depolarisation will occur 
even across a single channel at 2.4 GHz. 

Prior to periastron, T — 46 was the only occasion on 
which we were able to derive an RM from the pulsar. At 8.4 
GHz, the value we derived is —15200 rad m~ 2 . As detailed 
in Section 2, the data are obtained in 20 min sections at the 
high frequencies followed by 20 min at the low frequencies. 
We therefore divided the 4.8 GHz data into four 20 min sec- 
tions and derived the RMs independently. Figure ^ shows 
the position angles obtained across the 4.8 GHz band dur- 
ing the third observation on T — 46 (Sep 1) and the fitted 
slope. The RMs measured on this date are extremely high, 
the highest astrophysical RM measurement as far as we are 
aware. There is also marginal evidence for a varying RM 
throughout the observing period. We presume that the RM 
increased even further after this date ensuring depolarisa- 
tion even at 8.4 GHz on Sep 15 (T - 32). 

On T + 18 there is a marginal detection of the pulsar 
but no RM or DM can be obtained. On T+ 21, only the two 
lowest frequencies were observed but polarisation is clearly 
seen at 2.4 GHz and we derive an RM of —5100 rad m~ 2 . 
On T + 24, the RM has dropped significantly, however, on 
T+29 the RM is large and positive but can only be detected 
in the first 20 min integration at 4.8 GHz as detailed in 
the table. We note that this is the only positive value of 
RM obtained. A similar result was obtained during the 1994 
periastron, where the RMs were consistently negative both 
before and after periastron. Finally, as late as T + 63 the 
pulsar appears depolarised at all but 8.4 GHz. During the 
first 20 min observation at 4.8 GHz we measured an RM of 
-7700 rad m -2 . 

There is only marginal evidence for a variable DM with 
epoch, with the errors quoted in the table being much larger 
than the typical 0.2 - 0.4 c m~ 3 pc obtained using pulsar tim- 



ing (Johnston et al. 1996). The mean of the values quoted 
is 146.8 cm^pc, almost exactly the result obtained from 
timing. We note that these results are consistent with the 
values obtained during the 1994 and 1997 periastra. In 1994, 
the ADM (i.e. the DM with respect to 146.8 cm~ 3 pc) was 
3.7 cm~ 3 pc by T — 28 and no cha nge in DM was seen post- 
periastron (lohnston et al. 1996). In 1997, the ADM was 
1.5 cm~ 3 pc on T — 48 and 4 cm~ 3 pc on T — 42. An up- 
per limit of 0.2 cm~ 3 pc to the ADM was obtained post- 



periastron (lohnston et al. 2001) 



From the current data we therefore set an upper limit 
to ADM on T - 46 of 3.0 cm" 3 pc and on T + 29 of 5.8 
cm _3 pc (i.e. twice the error bars), although from the above 
discussion the ADM is likely to be much smaller on the 
latter date. The magnetic field parallel to the line of sight, 
By in /jG, can be obtained with knowledge of the RM and 
ADM by 



1.232 



RM 
ADM 



(•») 



where RM and ADM are in conventional units. The lower 
limits of B\\ on T - 46 (Sep 1) and T + 29 (Nov 15) are 
therefore 6.2 mG and 1.6 mG respectively. 
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Date (+T) Time Rotation Measure (rad m~ 2 ) DM 

UT 1384 MHz 2496 MHz 4800 MHz 8640 MHz cm" 3 pc 



Sep 01 


-46.4 


00:45 


- 01:05 






-14800±1800 








01:35 


- 01:55 






-6900±5000 








02:25 


- 02:45 






-12500±1400 








00:45 


- 03:45 








-15200±3500 148.7±1.5 


Nov 07 


21.3 








-5100±360 






Nov 10 


24.4 






-450±60 


-460±180 




145.7±2.9 


Nov 15 


29.5 


22:30 


- 22:50 






+7700±2100 


148.9±2.9 


Nov 20 


34.3 






-360±5 


-350±120 




144.1±3.7 


Nov 26 


40.3 






-675±30 


-750±100 




144.4±2.3 


Dec 11 


55.3 






-101±9 


-75±30 




148.4±1.6 


Dec 19 


63.4 


21:05 


- 21:30 






-7700±2400 


146.9±2.4 


Dec 26 


70.4 


19:25 


- 21:00 


-94±14 












21:00 


- 23:30 


-50±7 


-102±50 




147.8±2.0 


Jan 08 


83.1 






-49±9 


-24±35 




147.6±1.6 


Feb 07 


113.3 






13±11 


130±80 




145.1±2.5 



Table 3. RM and DM of the pulsar at various epochs averaged over the two pulse components. The errors are 1 a. Where no time range 
is given the entire observation was used. 




4 DISCUSSION 

A comparison of the light curves of 1997 and 2000 appears 
in Figure The same basic behaviour is seen for both sets of 
observations. There is a clear peak before and after perias- 
tron with a plateau in between. The transient emission also 
persists until late times in both. However, the initial 'blip' 
in the light curve seen at T — 21 in 1997 was not repeated 
in 2000. Johnston et al. (2001) identified this as a tran- 
sient event possibly associated with the pulsar 'splashing' 
into the disc. Also, in 2000, the first peak occurs earlier and 
the second peak later than in 1997. This is consistent with 
the longer lasting eclipse of the pulsed emission. The other 
striking feature is that the second peak is significantly lower 
in 2000 than in 1997. In broad terms, therefore, the model 
proposed after the 1997 periastron is repeated again. The 
pulsar enters the disc of the Be star near T — 20, electrons 



are accelerated in the interface between the pulsar wind and 
the disc material. This occurs over the few days that the 
pulsar takes to cross the disc. When the pulsar leaves the 
disc, the wind bubble remains behind and decays through 
synchrotron losses. During the second crossing of the disc 
near T + 20 the same process happens again causing the 
second peak seen in the light curve. 

In Fig |^ we again show the light curve at 1.4 GHz. The 
decay of the second peak is clearly not linear, but follows 
a more power-law like behaviour. This is the expected form 
of the flux density for adiabatic expansion of a synchrotron 
bubble. We can write the synchrotron flux density, S v , from 
an expanding bubble as 

S v cx L 2a ~ 2 S 1 "" v a (6) 
where L is the radius of the bubble, B is the magnetic field, 



© 2002 RAS, MNRAS 000, [j^§ 



The 2000 Periastron Passage of PSR B 1259-63 7 




-40 -20 20 40 60 80 100 120 
Time from periastron (days) 

Figure 6. Flux density of the impulsed emission at 1.4 GHz. 
Superposed is the solution to the model described in the text. The 
model includes a rise time of 10 days and adiabatic expansion for 
each of the synchrotron bubbles created in the disc crossings. The 
dashed lines indicate the light curve for each bubble individually. 
The solid line is the sum of both light curves. 



and v and a are the frequency and spectral index respec- 
tively. Imagine the synchrotron bubble is formed at an ini- 
tial distance Ro from the Be star and with an initial radius 
Lo- We then assume that the bubble moves outwards in the 
disc of the Be star with a velocity V. Hence 



R(t) = R + V(t- to) 



(7) 



where to is the start time of the adiabatic phase. The change 
in R is evidently R(t)/Ro. We assume that the density in 
the disc follows a power-law with distance from the star i.e. 
p(R) = CR~ n . It then follows that 



R + V{t-to) 
Ro 



(8) 



As the size of the bubble is dictated by the ram pressure 
between the pulsar wind and the disc material then L oc 
p~ 1//2 and so 



L_ = f Ro + V(t-t ) 
Lo [ Ro 



n/2 



(9) 



The magne tic field in the disc is radial at the d istances rel- 
evant here (Melatos, Johnston & Melrose 1995) and so 



B_ 



Ro + V{t-tp) 
Ro 



(10) 



Substituting these last two equations into equation q above 
we finally obtain 



S v = So 



Ro + V(t-t ) 
Ro 



(a-l)(l+n) 



(11) 



We model the data as follows. We set a — —0.6 from 
the observations and assume n — 3 (Waters et al. 1991). We 



assume that to for the first encounter is T — 8 and for the 
second is T + 21. The distance from the Be star at these 
epochs are 34 R, and 42 R* respectively. The interaction 
between the pulsar wind and the Be star's disc continues for 
~10 days. During this time electrons are injected continu- 
ously and a linear rise in the flux density provides a good fit 
to the data. We then fit the entire data set using the func- 



tional form of equation [ll] and with So for both encounters 
and V as free parameters. The solution yields So of 40 and 
24 mjy respectively (at 1.4 GHz) and V — 12 km s" 1 and 
this solution is shown in Fig ^| The ratio of the peak fluxes 
is in agreement with the relative distance from the Be star 
of the two start epochs. We note that although the fit is 
generally good it falls well below the measured flux density 
of the last data point at T+ 113. It is also evident that some 
unpulsed emission is present as early as T — 46 at least at 
low frequencies. There therefore appears to be some low- 
level emission unrelated to the disc crossings. This emission 
could arise with the interaction between the pulsar wind 
and dense clumps in the (more tenuous) Be star wind. We 
repeated the fitting process but in addition we had a con- 
stant offset of 4 mjy throughout the entire observing period 
to represent the low-level emission. This naturally provides 
an improved fit to both very early and very late times. How- 
ever, V is largely unchanged at 15 km s _1 . 

There is a strong similarity between the system under 
discussion here and the LSI+61°303 binary. There, an un- 
seen neutron star orbits a Be star in a highly eccentric orbit 
with a period of only 26 days. Radio outbursts show a steep 
rise followed by a power-law like decay. Paredes et al. (1991), 
in a more complete treatment, modelled the radio light curve 
of LSI+61°303 in a similar fashion to undertaken here - after 
a prolonged period of particle injection, they allowed the re- 
sultant synchrotron bubble to expand adiabatically. In their 
case, the light curve decays after only ~10 days, and they 
derive an expansion velocity of 400 km s _1 . In our case, the 
light curve decays on a much longer timescale which natu- 
rally implies a much lower expansion velocity. A velocity of 
only 12 km s _1 is surprisingly low for the outflow of disk 
material at this distance from the Be star. Waters et al. 
(1987) show that the velocity is likely to be ~100 km s _1 
and indeed our optical observations of this star indicate a 
velocity of ~80 km s _1 at 20 It is not clear, therefore, 
why our derived velocity is so low. 

Synchrotron loss times, set by the magnetic field in the 
bubble, are also relevant but on a longer timescale than 
the adiabatic losses. As the synchrotron bubble evolves, the 
break frequency in the emission spectrum moves to lower 
frequencies. Once the break frequency becomes lower than 
the observing frequency, the flux density at that frequency 
should decline rapidly. The time, tb, in days, at which the 
break frequency equals the observing frequency can be ap- 
proximated as 



t b ~ 315 B 



-3/2^-1/2 



(12) 



where B is the magnetic field strength in Gauss and v is in 
GHz. Given that no further injection of electrons occurs af- 
ter T + 22 and that the emission at 4.8 GHz appears to have 
decayed by T + 90 we can therefore estimate the decay of 
the emission at 8.4, 2.5, 1.4 and 0.84 GHz should occur 73, 
118, 150 and 185 days after periastron respectively. This is 
consistent with the lack of emission at 8.4 GHz from T + 63 
onwards and the detection of unpulsed emission at 1.4 and 
0.8 GHz at least as late as T+ 113 and T+ 132. Using Equa- 
tion [l2| the magnetic field in the disc where the synchrotron 
bubble is created is ~1.6 G. This is consistent with the result 
obtained by Ball et al. (1999). Note that, as expected, this 
value of the magnetic field in the disc significantly higher 
than the value we measured for the wind of the Be star. 
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5 CONCLUSIONS 

A series of observations of the 2000 periastron of the 
PSR B1259-63/SS2883 system were made at the ATCA at 
four frequencies. The pulsar binning mode available with the 
ATCA correlator meant that pulsed and unpulsed emission 
could be separated and we therefore obtained light curves 
for the pulsar and unpulsed emission independently. We were 
also able to obtain RMs and DMs for the pulsar. 

The main features in the light curve from the 2000 pe- 
riastron are similar to the 1997 periastron, providing confir- 
mation of our earlier model. It does appear, however, that 
the pulsar eclipse lasted longer and the unpulsed emission 
began and finished later in 2000 than in 1997. This indicates 
that the disc has become larger and/or denser in 3.5 yr. As 
we obtained more data at late times than in the previous 
periastron encounters, we found that adiabatic decay of the 
synchrotron bubble provided the best fit to the data. The 
expansion rate of the bubble, ~12 km s _1 , is significantly 
lower than expected but the magnetic field value of 1.6 G 
consistent with measurements in other Be stars. 

In addition, at T — 46 the pulsar had the highest mea- 
sured astrophysical RM of — 14800 ± 1800 rad m~ 2 implying 
a magnetic field of at least 6 mG in the Be star wind. 
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